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Abstract 

By employing a newly developed nonadiabatic dynamical simulation method, which is a com- 
bination of classical molecular dynamics (MD) and the adaptive time-dependent density matrix 
renormalization group (TDDMRG), we investigate the dynamics of charge carrier transitions in 
irans-polyacetylene (PA) with the inclusion of both electron-phonon and electron-electron inter- 
actions. The calculations are performed within a modified Su-Schrieffer-Heeger (SSH) model in 
which electron-electron interactions are taken into account via the combination with extended Hub- 
bard model (EHM). We find that removing an electron from a trans-VK chain bearing a positively 
charged polaron leads to the formation of a pair of charged solitons. Furthermore, we study the 
effect of electron-electron interactions on such charge carrier transitions in trans-FK. Our results 
show that increasing the on-site Coulomb interaction U and the nearest-neighbor Coulomb repul- 
sion V will not change the qualitative behavior of the transition from a polaron to a soliton pair 
in the evolution process but will quantitatively reduce the moving velocities of the both formed 
solitons significantly and change the conditions for the soliton collisions. 
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I. INTRODUCTION 



Since it was discovered in the 1970s that electrical conductivity of irans-polyacetylene 
(PA) can be improved significantly through charge injections or photoexcitations pj-y], con- 
jugated polymers have received sustained attentions from both academic and industrial 



ications, varying from 
4-9^ 



researchers because of their great numbers of conducting plastic app 
light-emitting diodes, field-effect transistors to photocells and lasers 

Different from the traditional conductors, the charge carriers in conjugated polymers are 
some self-localized nonlinear excitations, such as solitons, polarons, or bipolarons, which 
are inherent features of quasi-one-dimensional conducting polymer associated with various 
structural and charge distortions, depending on the properties of the polymer chain as well as 
on the concentration of charge doping. j8-12| In trans-FA, the most important quasi-particles 
responsible for the charge transport are solitons serving as domain walls to distinguish the 
two opposite bond length alternation patterns which might coexist in undisturbed trans- 
PA because of the two-fold ground state energetical degeneracy. A soliton may either be 
neutral with spin ±1/2 or present charge Q = ±e without spin. For the other conjugated 
polymers without ground state energetical degeneracy, polarons and bipolarons are the basic 
charge carriers for the charge transport in these materials. A polaron may present charge 
Q = ±e with spin ±1/2 and a bipolaron presents charge Q = ±2e with no spin. The 
studies of the transitions among these different charge carriers are therefore of great interests 
for the purpose of modulating electronic preperties of conjugated polymers under different 
circumstances. 

Recently, dynamics of the transitions between the different charge carriers in various 
conjugated polymers have been theoretically simulated by Silva et a/.jisl [l^ It was sug- 
gested that in conjugated polymers without ground state energetical degeneracy the most 
energetically favorable transition is a direct single-polaron to bipolaron transition while in 
trans-PA the most energetically favorable transition is that to a pair of charged solitons. 
However, in Silva et aPs dynamical simulations of charge carrier transitions electron-electron 
interactions are either completely ignored by only using Su-Schrieffer-Heeger (SSH) model 
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12| or merely partially included by adopting the Hartree-Fock (HF) approximation. 



Actually, many theoretical works have found that electron-electron interactions and elec- 
tron correlations play a crucial role in determining the nature of electronic excitations in 
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conjugated polymers. 15N35| For example, electron correlation is essential in determing the 
correct order of the lowest-lying vr — )■ vr* excitations in trans-PA. jsS] Conjugated polymers 
behave as quasi-one dimensional systems owing to their strong intramolecular interactions 
and rather weak intermolecular interactions and electron-electron interactions are accord- 
ingly weakly screened. However, accurate simulations of large quantum many-body systems 
by conventional post-HF methods are still too expensive due to the exponential growth of 
the number of freedom degrees. Fortunately, recent theoretical studies on the charge carrier 



transportation dynamics in conjugated polymers by our group and Wu et a/'s group [36|-[39| 



have shown that the combination of classical molecular dynamics (MP) and the adaptive 



40|-|42| provides a new 



time-dependent density matrix renormalization group (TDDMRG) 
efficient and accurate approach to perform real time nonadiabatic dynamical simulations 
of large quantum many-body systems including both electron-lattice and electron-electron 
interactions. In this paper, we apply such a method to simulate the charge carrier tran- 
sition dynamics in trans-PA under a photo-ionization of an electron to study the effect of 
electron-electron interactions on such dynamic processes. During the simulation, the time- 
independent and time-dependent Schrodinger equations are solved within the framework of 
SSH model modified to include electron-electron interactions via a combination with the 
extended Hubbard model (EHM). 

The paper is organized as follows. In Sec. II, we present the SSH model modified to 
include electron-electron interactions via EHM and describe the numerical method briefly. 
In Sec. HI the dynamics of the charge carrier transition dynamics in trans-PA under a 
photo-ionization of an electron will be discussed. A summary is given in Sec. IV. 



II. METHODOLOGY 
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12| with a combination with 



We use the well-known and widely used SSH Hamiltonian 
the extended Hubbard model (EHM) to describe both the vr-electron part and the lattice 
backbone with the inclusion of both the electron-lattice interactions and electron-electron 
interactions: 

H{t) = H,i + Hiatt. (1) 



The TT-electron part includes both the electron-phonon and the electron-electron interac- 



4 



tions, 



Hel = - ^ tn,n'{c^>^^Cn,a + h.C 
<n,n'>,cr 

u 

^1 



(2) 



1 1 



<n,n'>,(T,(T' 



where tn,n' is the hopping integral between the n-th site and the n'-th site, while U is the on- 
site Coulomb interaction and V denotes the nearest-neighbor electron-electron interaction. 
< n,n' > denotes the summation is over nearest-neighboring sites. 

Because the distortions of the lattice backbone are always within a certain limited extent, 
one can adopt a linear relationship between the hopping integral and the lattice displace- 



ments as 



n,n+l 



to — (y{un+i — Un) |l0l-ll2|. where to is the hopping integral for zero 



displacement, u„ the lattice displacement of the nth site, and a is the electron-phonon 
coupling. 

Because the atoms move much slower than the electrons, we treat the lattice backbone 
classically with the Hamiltonian 



K 



latt 



M 



Ur. 



(3) 



where K is the elastic constant and M is the mass of a site, such as that of a CH monomer 
for trans-PA. 

The model parameters are those generally chosen for trans-^^ 



A: to=2.5 eV, a=4.1 eV/A, 



K=21 eV/A2, M=1349.14 eVfsVA^ a=1.22 A, V = ?7/2.[l0|-ll3| We notice that there are 
some debates on the accuracy of some of the parameters. For example, Ehrenfreund, et al. 
proposed to use K=46 eV/A^ instead of 21 eV/A^ to fit the resonant Raman scattering 



experiment 43| , and an additional linear term was suggested to be included into to Eq. [3] 
for the more proper description of the elastic energy of the C-C a backbone j^l since a is 
assumed to be the equilibrium lattice spacing of the undimized chain, including both a- and 
TT-bonding. Fur the purpose of making contact to other theoretical studies on charge carrier 



transition problem 



13 



14l |. we choose the current parameter schemes. In order to prevent 



the contraction of the chain, which might result from our simplified lattice elastic energy 
form, we supplement a constraint of fixed total length of the chain. 
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The model chain is initially in its ground state configurations for the system containing a 
charged polaron. Therefore, we need firstly optimize the geometrical configurations of such a 
singly charged model chain. Due to the exponential growth of the number of degrees of free- 
dom in quantum many-body systems, the exact solution of our used SSH-EHM model in the 
complete Hilbert space is apparently not feasible for our system. Therefore, we perform ap- 
proximate solutions by virtue of the density matrix renormalization group (DMRG) method 
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45| . DMRG, firstly proposed by White in 1992 4^, |45|, which uses the eigenvalues of the 
subsystem's reduced density matrix as the decimation criterion of Hilbert space, has been 
shown as an extremely accurate technique in solving one-dimensional strongly correlated 



system with economic computational costs. 



42| 



Both strong electron-electron and electron-phonon integrations are main features of con- 
ducting polymers and the important reasons for the conjugated polymers to present novel 
photoelectronic properties. Therefore, it is obviously not reasonable for the backbone of 
conjugated polymers to be frozen in dynamic processes with a time of more than sev- 
eral femtoseconds. On the other way, if all the carbon-carbon bond lengths are fixed in a 
undimerized way, no solitons will be formed because there are no degenerate ground state 
possibilities. There might be also some charge density waves in such conditions, however, 
they will move forth and back much more quickly than the normal solitons and of course 
they are not the real case in conjugated polymers. 

For the purpose of performing real-time simulation of both the evolution of quantum vr- 
electron part and the classical movement of the chain backbone, we adopt a newly developed 
real-time simulation method in which classical molecular dynamics is combined with the 



adaptive TDDMRG |40H42| . The main idea of this method is to evolve the vr-electron 
part by the adaptive TDDMRG and move the backbone part by classical MD iteratively 
and nonadiabatically to include nearly all relevant electron-phonon and electron-electron 
interactions. Details about such the combined TDDMRG/MD method can be found in the 



recent papers 



36|-|39| 



III. RESULTS AND DISCUSSION 

The simulated system consists of a trans-PA chain with 100 (CH) monomers, where a 
positively charged polaron is initially localized in the middle of the chain. After the time 
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begins to evolve, one more electron is removed through photoionization. Then the dynamics 
of the charge carrier transitions in trans-PA within 400 femtoseconds (fs) is simulated by 
virtue of classical MD combined with the adaptive TDDMRG. 

Recent works have carefully studied the accuracy performance of the combined TD- 
DMRG /MD method and verified its reliability with suitable adopted parameters and it was 
found that one can achieve reliable and converged numerical results if the DMRG trunca- 
tion weight ep is taken small enough (e^ < 1.0 x 10~^) and suitable time step At values 



around 0.01~0.05 fs are chosen. 
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All our numerical results presented in this work are 
well-converged values calculated with Cp = 1.0 x 10^'' and At=0.05 fs. 

A. Structure of the initial state 

At t = 0, the singly charged model trans-PA chain with 100 (CH) monomers is located at 
its ground state. Its optimized geometric structure and charge distribution pattern at U=2.0 
eV and V^=1.0 eV are shown in Fig. [T]and Fig. |2j Considering there are strong bond length 
alternations and accordingly strong charge density oscillations due to significant electron- 
phonon and electron-electron interactions, for the purpose of presenting a better visualization 
of the results with more smooth curves, we use a mean charge density Pn(t), derived from the 
charge density Pn{t), and the staggered bond order parameter r^ft) by taking the average 



values of nearest-neighboring sites to analyze the simulations 



46| 



- /,^ _ Pn-l(t) + 2pn{t) + pn+l{t) / 
r4t) = (_i)n^n-l(t)+2Mt)+^n+l(t)^ 

where Pnit) = 1 — J^ai"^ i^)\i^n a(^n,a\'^ (t)) and Un is the lattice displacement of the nth 
site as mentioned earlier. One may find typical characteristics of polarons in conducting 



polymers [12| from the geometric structure and charge distribution pattern before photoion- 
ization as illustrated in Fig. [T] and Fig. [21 This implies that the initial singly charged model 
chain is bearing a polaron defect. The model chain is then assumed to undergo further 
vertical photoionization in the spirit of the Franck-Condon principle, which means one more 
electron is taken away from the original singly charged model system immediately at the 
unchanged ground state geometry of singly charged model chain. The removal of the sec- 
ond electron originally located at singly occupied molecular orbital (SOMO) makes the new 



doubly ionized state lying at higher energy level than the ground state of the original singly 
charged system. Such kind of energy difference is determined to be about 0.84 eV by our 
calculations with [7=2.0 eV and ^=1.0 eV, implying that such doubly ionized state locates 
between the second and third lowest doublet excited states of singly charged trans-PA, whose 
vertical excitation energies were predicted to be about 0.58 and 0.91 eV respectively with 



Pariser-Parr-Pople (PPP) calculations 301 ]. 



(Fig. [T] about here) 
(Fig. [2] about here) 



B. General picture of charge carrier transitions in trans-PA 

(Fig. [3] about here) 
(Fig. m about here) 

Now let us show the general time evolution pictures of charge density and site displace- 
ment in the charge carrier transition process in trans-PA. The time evolution pictures of 
the mean charge density Pn{t) and the staggered bond order parameter r„(t) are shown in 
Fig. [3] and Fig. |4] respectively. 

As we can see from Fig. [31 the removal of an electron from the trans-PA chain which 
initially bears a positively charged polaron induces larger charge densities in the central area 
of the chain. Then the induced large charged densities in the middle chain tend to disperse 
to the nearby areas toward the both ends of the chain while the geometrical distortions 
will also evolve as illustrated in Fig. |4] due to the electron-phonon interactions. Because 
the two opposite bond length alternation patterns are energetically degenerate in trans-PA, 
it is energetically more favorable for the new excitation to break into a pair of charged 
solitons than a pair of polarons. It could be found that a pair of charged solitons with two 
charge peaks have been formed after about 20 fs from the pictures of the charge density and 
geometrical structure in Fig. [3] and Fig. |H Clearly, the two charged solitons repel each other 
because they are with the same charge sign. From Fig. [3] we may find that the moving charge 
peaks do not really change at all (but for the positions) during the repulsion processes of the 
two new defects which shows that the transition from polaron to charged soliton pair has 
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been finished and the two newly formed charged solitons can move relatively freely after 20 
fs. From Fig. H] it is found that during the movement of the two charged solitons the center 
positions of the solitons move as nearly strait lines with time evolving which implies the 
solitons are moving with a nearly constant velocity. It could also be found that the charge 
distortion is always coupled with the geometrical distortions very well during the entire 
moving process of the soliton which verifies that the soliton structure is an inherent feature 
of trans-PA. In Fig. H] we can notice some wave grains appear after the solitons begin moving. 
That happens because there is always a long-lasting oscillatory "tail" appearing behind the 
soliton defect center during the soliton moving process. 36|] This "tail" is generated by the 
inertia of those monomers to fulfill energy and momentum conservation and it absorbs 
the additional energy, preventing the further increasement of the soliton velocity after a 
stationary value is reached. 

As illustrated in Fig. [3] and Fig. HI the two charged solitons seem to be able to move along 
the chain back and forth independently after they are formed at around 20 fs. They get 
reflections when they encounter the boundaries of the chain at around 140 fs, and then they 
return back and get the collisions when they meet each other at around 280 fs. Afterwards, 
the charged soliton pair will repeat the processes of repulsion, independent moving and 
collision. 



C. Effect of electron-electron interactions on the charge carrier transitions in 
trans-PA 

In order to study the influence of electron-electron interactions on the charge carrier 
transitions in trans-PA, we perform various simulations for the charge carrier transition 
dynamics with different adopted values for the electron-electron interaction strengths. Time 
evolution pictures of the mean charge density pn{t) and the staggered bond order parameter 
r„(t) simulated with different the on-site Coulomb interaction U values and the nearest- 
neighbor Coulomb repulsion V values are displayed in Fig. \5\ and Fig. El 

(Fig. O about here) 

(Fig. [6] about here) 
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From Fig. E] and Fig. E] we can clearly see that the evolution pictures for both the 
charge densities and the geometrical parameters with different electron-electron interaction 
strengths are qualitatively very similar. This means the varying of U and V will not change 
the qualitative behavior of the transition from a polaron to a soliton pair in trans-PA after 
photo-ionization: removing an electron from the chain initially bearing a positively charged 
polaron leads to the formation of a pair of charged solitons and then the formed two soli- 
tons will move back and forth along the chain independently as well as two solitons will get 
reflections at the chain boundaries and get collisions when encountering each other in the 
chain center region. 

However, we may found that increasing both U and V will quantitatively reduce the 
moving velocities of the two newly formed solitons as illustrated in Fig. |5] and Fig. [61 Such 
an observation is in good agreement with our previous studies on the soliton transport 
dynamics in trans-PA under an external electric field 36||, in which we have found that in 
most cases {U > 0.9 eV) increasing U will induce the reduction of electron delocalization 
and accordingly it will suppress the transport of the charged soliton. Our previous work 
has also shown that increasing U when U < 0.9 eV will favor the charged soliton transport 
while increasing V when V is small will be very unbeneficial to the charged soliton transport 
because it will induce a more localized defect distribution. 



Therefore, it is reasonable to 



see increasing both U and V at the same time will suppress the movement of two newly 
formed charged solitons, as illustrated in Fig. |5] and Fig. |6l 

At the same time, we may also notice the varying of U and V will change the conditions 
for the collision of the two newly formed charged solitons from Fig. [5] and Fig. [61 With 
the enhanced U and V, the minimal distance between the two charged solitons during the 
collision process will be significantly increased, from less than 10 units in the case with 
U=0.0 eV and V=0.0 eV to about 30 units in that with U=3.0 eV and ^=1.5 eV. This 
can be easily understood for the fact that the enhanced electron-electron repulsions will 
make the repulsion between two charged solitons stronger and accordingly the two charged 
solitons can hardly move furthermore to touch each other completely under the case with 
large U and V values. 

Actually, the case with U =0 and ^=0 in our presented results is the same as what has 
been studied in Silva et al work with pure SSH modeljl^. As we have disclosed before, 
the further consideration of the on-site Couloumb repulsions U and the nearest-neighbor 
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electron-electron interactions V will not change the qualitative behavior of the charge carrier 
transition in trans-PA after photo-ionization, but will reduce the moving velocities of the 
two newly formed solitons significantly as a result of decreasing electron delocalization. 

In this work, we use V values varying from to only 1.5 eV which might be not so realistic 
in view of the spectroscopic properties of ethylene. The further increase of nearest-neighbor 
interactions V will strengthen the electron-hole attractions and accordingly induce much 
more delocalized soliton picture according to previous studies js^. Considering our current 
model chain has only 100 sites and the increase of U and V do not change the qualitative 
behavior of the charge carrier transition in trans-PA at least within this study, we don't use 
larger U and V values in this work for the fact that it is difficult to present the pictures of 
two moving and much more delocalized solitons very clearly in a short chain. However, we 
are expecting to simulate larger model systems with more realistic electronic integrals with 
more powerful computational facilities in the future works. 

IV. SUMMARY 

For a model trans-PA chain initially holding a polaron defect, described with a modified 
SSH model in which electron-electron interactions are included via the combination with 
EHM, we have studied the dynamics of charge carrier transitions in such a trans-PA chain 
after photo-ionization. Nearly all relevant electron-phonon and electron-electron interactions 
have been fully taken into account by solving the time-dependent Schrodinger equation for 
the vr-electrons and the Newton's equation of motion for the backbone monomer displace- 
ments with the combination of the adaptive TDDMRG and classical MD. 

It is found that removing an electron from the trans-PA chain bearing a positively charged 
polaron leads to the formation of a pair of charged solitons and then the formed two soli- 
tons will move back and forth along the chain independently. These two solitons will get 
refiections at the chain boundaries and get collisions when encountering each other in the 
chain center region. 

Furthermore, we study the infiuence of electron-electron interactions on such charge car- 
rier transitions in trans-PA. Our results show that increasing the on-site Coulomb interaction 
U and the nearest-neighbor Coulomb repulsion V will not change the qualitative behavior of 
the transition from a polaron to a soliton pair in the evolution process. However, increasing 
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both U and V will quantitatively reduce the moving velocities of the two newly formed 
solitons significantly as a result of decreasing electron delocalization and increasing U and 
V will also change the conditions for the soliton collisions. 
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Figure Captions 

FIG. □ Geometric structure at the initial state. {U=2.0 eV, ^=1.0 eV) 

FIG. [2] Charge distribution picture before and after vertical photoionization at the 
initial state. (f/=2.0 eV, V=1.0 eV) 

FIG. [3] Time evolution of the mean charge density p„(t). (f/=2.0 eV, 1^=1.0 eV) 

FIG. m Time evolution of the staggered bond order parameter r„(t). (f/=2.0 eV, 
V=1.0 eV) 

FIG. E] Time evolution of the mean charge density Pn{t) with various electron-electron 
interaction parameters. 

FIG. [6] Time evolution of the staggered bond order parameter r„(t) with various 
electron-electron interaction parameters. 
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FIG. 1: Ma et al, J. Phys. Chem. 
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FIG. 3: Ma et al, J. Phys. Chem. 
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FIG. 5: Ma et al, J. Phys. Chem. 
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FIG. 6: Ma et al, J. Phys. Chem. 
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